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We examine nonclassical properties of the field states generated by applying 
the photon annihilation-then-creation operation (AC) and creation-then- 
annihilation operation (CA) to the thermal and coherent states. Effects of 
repeated applications of AC and of CA are also studied. We also discuss 
experimental schemes to realize AC and CA with a cavity system using 
atom-field interactions, (c) 2009 Optical Society of America 

OCIS codes: 270.5290, 270.5585. 
1. INTRODUCTION 

In many applications of quantum information processing that employ radiation fields, 
it is often desirable that the fields possess nonclassical nature such as entanglement 
and squeezing. From a theoretical point of view, perhaps the simplest way of generat- 
ing nonclassical field states is to apply the photon creation operation a+ to classical 
states such as the thermal and coherent states [1]. These nonclassical field states, 
referred to as the single-photon added coherent state [2] and the single-photon added 
thermal state [3], respectively, have been realized experimentally [4-6]. 

In this paper, we investigate nonclassical properties of the field states that re- 
sult when the photon annihilation-then-creation operation (AC) a^a and the photon 
creation-then-annihilation operation (CA) aa^ are applied to the thermal and coher- 
ent states. The noncommutativity of the photon operations, [a, a^] = 1, breaks the 
symmetry between AC and CA. It thus is of interest to analyze and compare the 
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properties of the states that result after AC and CA. In fact, a direct observation of 
the noncommutativity of the operators a and a"*" has been achieved by experimen- 
tally comparing such two states and proving that they differ from each other [7]. 
An experimental scheme to directly prove the commutation relation [a, a'^] = 1 has 
also been proposed recently [8]. Here we look in particular at the possibility and the 
degree of nonclassicality exhibited by the states produced by AC and CA. We also 
investigate properties of the states obtained when AC and CA are applied repeatedly 
to the thermal and coherent states. 

The paper is organized as follows. In Sec. II we briefly review the effect of the 
photon creation and annihilation operations a"*" and a. In Sec. Ill we study properties 
of the states produced by AC and CA operated on the thermal state and the coherent 
state. We also study the effect of repeated applications of AC and CA to the thermal 
and coherent states. Experimental schemes to realize AC and CA with a cavity system 
are described in Sec. IV. Finally, Sec. V presents a discussion of our results. 

2. CREATION AND ANNIHILATION OPERATIONS 

The creation (annihilation) operator a"*" (a) acting on a photon number state \n) 
results in an increase (a decrease) of the photon number by one. As the creation (an- 
nihilation) operator is not a Hermitian, however, it cannot be expected to always play 
the role of physically adding (subtracting) one photon. Let us consider an arbitrary 
pure field state \ip) = Xl^o ^nln), where J2'^=o \^nf = 1- The state that results after 
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application of the operator a"*" (a) to the state 1-0) is \tlj)c = ^{^'''l'^)} 



y/{^\aa+\'^) 



= N{a\'(Jj)} = — ^=M==), where signifies that normahzation is to be per- 



formed on the unnormahzed state inside the curly bracket. One easily finds 

{h)c - {n) 



{n) + 1 

in) A -in) = (2) 

where (n), {h)c and (n)^ refer to the mean photon number of the state \ip), \4')c 
and \iP)a, respectively, and (Ari)^ is the photon-number variance of the state I'i/;) 
[(An)^ = (n^) — (n)^]. Eqs. (1) and (2) indicate that, only when the state \ifj) has zero 
photon-number variance, the creation (annihilation) operation increases (decreases) 
the mean photon number by one. The mean photon number increases at least by 
one under the action of the creation operator and decreases at most by one under the 
action of the annihilation operator. In fact, the mean photon number can even increase 
upon application of the annihilation operator, which occurs when (An)^ > (n). This 
tendency of the mean photon number to increase beyond one's naive expectation has 
its origin in the relation a~^\n) = y/n + l\n+ 1) and a\n) = y/n\n — 1); a large photon 
number carries a greater weight in the state obtained after the a"*" or a operation than 
in the original state. 

The operations a"*" and a can have more profound effects upon the state they are 
applied to than simply changing the number of photons. In fact, the states that result 
after the operation a"^ and a often exhibit significantly different physical properties 



from the original state and from each other. It is well known that a creation operation 
on thermal and coherent states yields nonclassical states, while the states remain 
classical upon an annihilation operation. 

3. ANNIHILATION-THEN-CREATION AND CREATION-THEN- 
ANNIHILATION OPERATIONS 

In this section we consider the effect of the annihilation-then-creation operation (AC) 
a'^a and the creation-then-annihilation operation (CA) aa~^ on thermal and coherent 
states. We further study the properties of the states that are obtained when AC and 
CA, respectively, are applied repeatedly. We look in particular at the possibility and 
the degree of nonclassicality exhibited by these states. 

3. A. Initial thermal state 

The thermal state is characterized by the density matrix 



where n = ^nu,/lT_i ■ The mean photon number {n) and the photon- number variance 
{Ahf = {h^) - {hf for the state (3) are given respectively by (n) = n and (An)^ = 
+ n, yielding the Mandel Q factor [9] Q = 1 = n. Application of the 




oo 



(3) 
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operation AC and CA, respectively, to the state (3) yields 



Pac = JS \a apa a] — 



Tr{a'^apa^a] 

-, oo „ 2 

1 rl^n^ 

^ > 7 , (4) 

l + 2rr 1 + n 1 + n ' 



= N{a(i^paa'^} 



n=0 

aa'^paa'^ 



Tr{aa+ paa'^} 



(l + 2rr)(l + n) ^ 

The Mandel Q factor for the states (4) and (5) can be obtained through straightfor- 
ward calculations. We show in Fig. 1 the calculated Q factor as a function of n, the 
mean photon number of the initial thermal state (3). It can be seen that, when the 
AC operation is applied to the thermal state, the Q factor becomes negative and thus 
the photon number distribution becomes sub-Poissonian for a sufficiently small initial 
mean photon number, i.e., for n < 0.6. The sub-Poissonian photon statistics implies 
nonclassical properties of the state, although the reverse is not necessarily true. Under 
the CA operation, however, the photon statistics remains super-Poissonian regardless 
of n. The figure indicates also that the difference in the Mandel Q factor between the 
two states pAc and pcA is more significant when n is smaller. The difference in the 
degree of nonclassicality between the two states (4) and (5) is illustrated in Fig. 2, 
in which the corresponding Wigner distribution functions [10, 11] are plotted for the 
case n = 0.57. Only the Wigner distribution for the state (4) takes on negative values 
in the vicinity of the origin. 

We now consider the states produced by repeated applications of a'^a or aa+ to the 



initial thermal state (3). A straightforward algebra yields 



'Vl+n^ n=0 ^ ' 

p'c^ = jv{(aa+)V(<i<i+)n 

where k can be any positive integer and Lis{z) = Xlfcli f polylogarithm. The 

mean photon number {{n)^^ and {n)%A) the variance ([(An)^]^^; and [(An)^]^^) 
for the states (6) and (7), can also be obtained by straightforward calculations. We 
obtain 

{n)\c = ^l:"'-ft(^ ={n)'cA + l, (8) 

2-lk \fA£^\2ik 



(9) 



We thus see that the mean photon number of state (6) is always larger by exactly one 
than that of state (7) (Note, on the other hand, the operator identity a~^a = aa^ — 1-) 
regardless of k. We also see that the photon number variance is exactly the same for 
the two states (6) and (7), regardless of k. 

In general, as the number k of operations a"*" a or aa^ is increased from 1, the mean 
photon number of the state increases rapidly. This is due to the fact already mentioned 
in Sec. II that a larger photon number carries a greater weight in the state after the a"*" 



or a operation than in the original state. It then is clear from Eqs. (8) and (9) that, for 
large k, the two states and possess almost identical photon-number statistics. 
In Fig. 3 we plot the Mandel Q factor as a function of n for the two states p'Xc^ 
and PcA^- The two curves are seen to be very close to each other. Close inspection 
reveals, however, that the curve for the state p'}^^ lies slightly below that for the state 
PcA^°5 indicating that the state p^c^ shows slightly stronger nonclassicality than the 
state p^A^ ■ For both states, the Mandel Q factor is negative and the photon-number 
distribution is sub-Poissonian if n < 0.57. That the two states p\(j and approach 
each other as k moves toward a large number can clearly be seen in Fig. 4, in which 
the fidelity between the two states, F = {Tr[^^/ p\cPcaV PacV^"^}"^ [l^] is plotted as 
a function of k, for the case n = 0.57. 

3.B. Initial coherent state 

The coherent state \a) has a Poissonian photon-number distribution, and thus Q = 
for the coherent state. Upon application of AC and CA, respectively, the coherent 
state \a) is transformed to the states 

|«)ac = N{a+a\a)} = ^=^^==, (10) 

{a\a+aa+a\a) 

\a)cA = N{aa^\a)} = ^=^^M==. (11) 

It should be mentioned that, since the coherent state is an eigenstate of the annihi- 
lation operator a, \a)Ac is given simply by \a)Ac = N{a'^\a)} and can be identified 
as the single-photon added coherent state. 



Fig. 5 shows the Mandel Q factor for the states (10) and (11) as a function of 
|a| , the mean photon number of the initial coherent state. It is seen that both states 
(10) and (11) have negative Q values regardless of As before for the case of the 
initial thermal state, the AC operation produces stronger nonclassicality than the CA 
operation, with the difference more significant when \a\ is smaller. In Fig. 6 we show 
the Wigner distributions for the states (10) and(ll) for the case \af = 0.57. One 
sees that both AC and CA operations lead to negative Wigner distributions, but the 
degree of negativity created by the AC operation is stronger. 

Let us consider the states produced by repeated applications of AC and CA, re- 
spectively, to the coherent state, given by 

= N{{a^a)'\a)}, (12) 
\a)'cA = N{{aa+)'\a)}. (13) 

Shown in Fig. 7 is the Mandel Q factor for the state \o()ac^ ^^"^ \'^)c^^ ^ function 
of I a I . One sees here that both states exhibit negative Q factors, but the state 
\^)ac^ shows a slightly stronger sub-Poissonian photon distribution than \(y)'^^. 
Comparing Fig. 5 and Fig. 7, one also sees that, as k is increased, the two states 
and approach each other. In Fig. 8 we plot the fidelity between the two 

states, F = \a_c{'^\^)ca\ 5 ^is a function of k, for the case |a| = 0.57. The fidelity 
approaches 1 as A; is increased. 
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4. PHYSICAL IMPLEMENTATIONS 

The photon creation operation a"*" on the coherent and thermal states was reahzed 
experimentally by Zavatta et al. using an optical system based on parametric down- 
conversion in a nonlinear crystal followed by conditional detection of a single photon 
in the idler mode [4-6] . The same group later achieved an experimental realization of 
the photon annihilation-then-creation operation (AC) a'^a and the photon creation- 
then-annihilation operation (CA) aa"*" on the thermal state using again an optical 
system, which provided direct observation of the noncommutativity of the operators 
a+ and a [7]. Here we briefly desribe an experimental scheme to realize AC and CA 
using a cavity system along the line suggested by Sun et al. [13] . 

In order to realize AC, we pass two atoms, atom 1 prepared in the lower state 
\g) and atom 2 in the upper state |e), in order through a cavity. Let us assume an 
arbitrary pure state = Yl'^=o^n\n) for the initial cavity field. (For simplicity we 
assume a pure state. Essentially the same analysis can be applied to obtain the same 
result even if the initial cavity field is in a mixed state.) A simple algebra yields that 
the state of the cavity field after the interaction of the cavity with the two atoms, 
conditioned upon the observation that atom 1 is found in state |e) and atom 2 in 
state l^f) after the interaction, is given by 

oo 

= ^iYl ^HV^gti) sin(v^(/t2)C„|n)} (14) 

n=l 

where g is the atom-field coupling constant, and ti and t2 are, respectively, the inter- 
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action time between atom 1 and the cavity and between atom 2 and the cavity. The 
state (14) can approximately be identified as the state 

oo 

I V)ac = N{a+a\^),} = iV{^ nCn\n)} (15) 

n=l 

if the interaction times ti and t2 are sufficiently short that {n)gti ^ 1 and 
\/Jn}gt2 <^ 1, where (h) is the mean photon number of the initial field state 
Repeated operations of AC, (a'^a)'', can then be realized by letting k pairs of atoms, 
each pair prepared in l^f) and |e), interact with a cavity, and postselecting the case 
where each pair is found in |e) and \g), subjected to the condition that all interaction 
times are sufficiently short. 

Similarly, the operation CA can be realized by letting two atoms, atom 1 prepared 
in |e) and atom 2 in \g), interact with a cavity field, conditioned upon atom 1 being 
found in \g) and atom 2 in |e) after the interaction and the interaction times ti and ^2 
being sufficiently short. In this case, the state of the cavity field after the interaction 
is 

oo 

IV')2 = N{J2 sm{Vn + lgh) sin(Vn + l^t2)C„|n)} (16) 

n=0 

which, in the short time limit a/ (n) + Igti ^ 1, a/ (n) + lgt2 ^ 1, approaches the 
state 

oo 

\ilj)cA = N{aa+\ilj),} = N{J2{n + l)C„|n)}. (17) 

n=0 

Repeated operations of CA, (da"'')'', can also be realized by using k pairs of atoms. 

11 



The AC and CA operations can be realized with the cavity scheme described above 
with high fidehties Fi = 11(^1-0)^01^ — 1 and F2 = |2('0|V^)ca|^ — 1, if the interaction 
times are sufficiently short. The problem, however, is that, in the short-time limit, 
the success probability of the scheme is low. Let Pi (P2) be the success probability of 
the scheme to realize AC (CA), i.e., let Pi (P2) be the probability to find atom 1 in 

|e) {\g)) and atom 2 in \g) (|e)) after the interaction. One can easily find that 

00 

Pi = ^sin2(v^^ti)sin'(v^(7t2)|C„|', (18) 

n=l 

00 

P2 = ^sin2(v/^rri^ti)sin2(v/^rTI^t2) (19) 

n=0 

In the short-time limit, Pi,P2 oc g^t\t\ <^ 1. For repeated operations {o^aY and 
{aa^Y , the success probabilities are extremely small in the short-time limit as Pi, P2 oc 

„4fcj.2fc-t2A; 
y ''1 ''2 ■ 

Instead of taking the short interaction times, one can choose ti and t2 in such a 
way that \J {n)gti = a/ {n)gt2 = | for AC and a/ (n) + Igti = a/ {n) + lgt2 = | 
for CA, where (n) is the mean photon number of the initial cavity field. This choice 
ensures that the desired state transitions \g) ^ \e) for both atoms 1 and 2 occur 
with high probabilities. When the atom makes a transition as desired, the cavity field 
gains or loses a photon as desired. One can thus hope that, with this choice of the 
interaction times, both the success probability and the fidelity may stay reasonably 
close to unity. Clearly, this strategy would work well if the initial cavity field has a 
sharp photon-number distribution centered at (n), which is the case for the initial 
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coherent state but not for the initial thermal state. In Fig. 9 we plot the success 
probability Pi and the fidelity Fi we computed as a function of |a| for the case when 
the initial cavity field state is a coherent state, where ti and ^2 are chosen to satisfy 
the relation a/ {n)gti = ^/(n}gt2 = f . One can see clearly that both the success 
probability and the fidelity stay close to unity for {n) > 50. The success probability 
P2 and the fidelity F2 show similar behavior, the only significant difference being that, 
as (n) ^ 0, Pi ^ as can be seen from Fig 9 but P2 ^ 1. The result shown in Fig.9 
suggests that, with this choice of the interaction times, even repeated operations of AC 
and CA can be realized experimentally with reasonably high success probability and 
fidelity, if the initial cavity field has a sufficiently sharp photon-number distribution. 

To assess the feasibility of the suggested experiment, we take, as an example, g = 
10^ Hz and (h) = 100. The interaction times ti and ^2 are then required to be ^1 = ^2 = 
—^^==- = 10~^sec. Since photon addition and subtraction should occur within these 
interaction times, we require the spontaneous emission decay rate 7 and the cavity 
field decay rate k be small compared with ^ = ^5 we require 7 <C 10^ Hz and 
K ^ lO^Hz, an experimental condition that can be met without too much difficulty. 

Photon annihilation-then-creation and creation-then-annihilation operations can be 
implemented experimentally using cavity-field atom interactions as described above 
as well as using optical means along the line of Refs. [4-6] . The cavity-QED method 
described here offers the flexibility of choosing interaction times in such a way that the 
desired states \iP)ac and \ip)cA can be generated with both high success probability 
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and high fidehty. The noncommutativity of the photon creation and annihilation 
operations can thus be verified with high success probabihty not easily achievable in 
the optical approach, if the cavity-QED method is adopted. 

5. DISCUSSION 

We have studied effects of the photon annihilation-then-creation operation (AC) a~^a 
and the photon creation-then-annihilation operation (CA) aa"*", and their repeated 
applications on the thermal and coherent states. The common feature of AC and CA 
is that they both work to increase the mean photon number unless the initial state 
is the number state. That this is the case can be understood by noting that, due to 
the relation a~^a\n) = n\n) and aa'^\n) = {n + a larger photon number carries a 

greater weight in the state obtained after the operation a^a or aa^ is applied than in 
the initial state. Details of exactly how the weight of each photon number is changed, 
however, are different for the two operations. Close inspection of this difference shed 
light on one of the main observations of this work, namely the observation that AC 
in general produces stronger nonclassicality than CA. It should be noted, in partic- 
ular, that the vacuum state, which is a part of the initial thermal or coherent state, 
disappears upon application of the a~^a operation, whereas it still remains after the 
aa+ operation. As a result, the state produced by the a+a operation has a signifi- 
cantly different photon number distribution from the initial "classical" distribution, 
especially if the vacuum state has a large weight in the initial state, i.e., if the ini- 
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tial thermal or coherent state has a small mean photon number. This explains why 
nonclassical properties are exhibited more strongly after AC than after CA and why 
the difference is more pronounced when the mean photon number of the initial state 
is small. If the operation a~^d or da~^ is applied repeatedly, the mean photon number 
of the state increases and accordingly the weight of the vacuum state becomes small. 
In such a case, the effect of the d'^a operation will not be much different from the 
effect of the dd'^ operation. This explains why the two states produced by applying 
the a^a operation k times and the ad"*" operation k times approach each other as k 
is increased. 

In this work, we have used the Mandel Q factor and the Wigner distribution as an 
indicator of nonclassicality. In a recent work [14], nonclassicality of the field states 
generated by photon creation-then-annihilation operations were investigated using 
the nonclassical depth [15] as a criterion for nonclassicality. We further note that 
our investigation here focuses on nonclassical effects related to sub-Poissonian pho- 
ton statistics. As photon-added thermal and coherent states have been shown to ex- 
hibit higher-order nonclassical effects beyond negativity of the Mandel Q factor and 
the Wigner distribution function [16], the states generated by photon creation-then- 
annihilation operations and by photon annihilation-then-creation operations should 
also exhibit such higher-order nonclassical effects. 
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FIGURE CAPTIONS 

Fig. 1: Mandel Q factor of the states pac (dashed curve) and pcA (sohd curve) vs. n, 
the mean photon number of the initial thermal state. 

Fig. 2: Wigner distribution of (a) the state pac and (b) the state pcA for the case 
n = 0.57. X = |(a + a*) and y = ^(a — a*). 

Fig. 3: Mandel Q factor of the states p^c^ (dashed curve) and p^'^ (solid curve) vs. 
n. 

Fig. 4: Fidelity between the states p\^ and vs. k for the case n = 0.57. 

Fig. 5: Mandel Q factor of the states |a)yic (dashed curve) and \a)cA (solid curve) 

vs. I a I , the mean photon number of the initial coherent state. 

Fig. 6: Wigner distribution of (a) the state \a)Ac and (b) the state \a)cA for the case 
\af = 0.57. X = |(a + a*) and y = ^{a — a*). 

Fig. 7: Mandel Q factor of the states \<y)'Ac^ (dashed curve) and Ict)^^^^ (solid curve) 
VS. \a\ . 

Fig. 8: Fidelity between the states \a)\Q< and la)^^ ^s. k for the case \af = 0.57. 
Fig. 9: The success probability Pi (solid curve) and the fidelity Fi (dashed curve) vs. 
|a| , the mean photon number of the initial coherent state. The interaction times ti 
and ^2 are chosen such that \a \ gti = \a \ gt2 = |. 
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